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ABSTRACT 

Rapidly rotating Neutron Stars in Low Mass X-ray Binaries (LMXBs) may be an in- 
teresting source of Gravitational Waves (GWs). In particular, several modes of stellar 
oscillation may be driven unstable by GW emission, and this can lead to a detectable 
signal. Here we illustrate how current X-ray and ultra-violet (UV) observations can 
constrain the physics of the r-mode instability. We show that the core temperatures 
inferred from the data would place many systems well inside the unstable region pre- 
dicted by standard physical models. However, this is at odds with theoretical expec- 
tations. We discuss different mechanisms that could be at work in the stellar interior, 
and we show how they can modify the instability window and make it consistent with 
the inferred temperatures. 
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1 INTRODUCTION 

Low Mass X-ray Binaries (LMXBs) were suggested as 
interesting sources o f Gravitational Waves (GWs) more 
than thirty years ago (|Papaloizou fc Pringldll978l : IWagonerl 
1984). In these systems, a compact object, which in the case 
of interest is a neutron star (NS), accretes mass from a less 
evolved low mass companion. The mass donor fills its Roche 
lobe, and matter is stripped from the outer layers and forms 
an accretion disc. The disc matter gradually loses angular 
momentum and spirals in, until it is eventually accreted 
by the NS. This process leads to angular momentum being 
transferred to the NS which can then be spun up to mil- 
lisecond periods in what is know as the "recycling" sc enario 
l|Alpar et al.lll982l : iRadhakrishnan fc Srinivasanll 19821 ) 

The main reason for invoking GW emission from these 
systems is the fact that the distribution of spin rates of 
both LMXBs and millisecond radio pu lsars (MSRPs) ap- 
pears to have a cutoff at around 730 Hz (|Chakrabartv et al.l 
2003) , which is w e ll below the centrifu gal break up limit 
I Cook et al.l 1 1994 lHaensel et all ll999l V This observation 
still holds true today, even as more systems have been 
added to the sample jPatrund [ioioh . Thus it is natural 
to seek a physical mechanism that can prevent NSs from 
spinning up further. The most obvious candidate is the ac- 
cretion process itself, as the interaction between the accre- 
tion disc and the star can lead to spin equilibrium if the 
system approaches a propeller phase and further accretion 
is centrifugally inhibited. This mechanism dictates a cor- 



relat ion between the magn etic field strength and accretion 
rate (|White fc Z hang 1993) i a problem which led to the pro- 
posal of several GW emission mechanisms that could gen- 
erate a s trong enough torque to set the spin equilibrium of 
LMXBs (fB ildstcn 1998; Andcrsson, Kokkotas & Stcrg ioulasl 
1 19991 : lUshomirskv. Cutler fc Bildstenl l20od : ICutlerl |2002| ). 
Although several authors have reassessed this problem 
jAndersson et a l. 2005; Ho. Macca rone fc Anderssonl boill : 
iPatruno. Haskell fc D'An gclo 2011), the question remains 
unresolved and current GW se arches are not sen sitive 
enough to give strong constraints (|Abbott et al.ll2010r i. 

The main GW emission mechanisms that could be at 
work i n accreting sy s tems are "mountains", either on the 



crust jBildstenlll998l: rUshomirsky. Cutler fc BildstenH2000|: 



Haskell. Jones fc Andersson 20061) or i n the core j Owenl 
19951 : 1 Andersson. Haskell fc Comerll2010|). deformations due 
to th e magnetic field of the st ar ~(|Cutlerll200"2l : lHaskell et all 
120081 : iMelatos fc Pavnell2005l ). and modes of oscillation of 
the st ar being driven unstable and growing to l arge ampli- 
tudes (|Andersson. Kokkotas fc Stergioulaslll999l ). 

We shall focus on the last, specifically the r-mode 
instability. An r-mode is a toroidal mode of oscillation for 
which the restoring force is the Coriolis force. It is particu- 
larly interesting be cause it is not on l y generically unstable 
to G W emission (|Anderssonl 1 19981 : iFriedman fc Morsinkl 
1998) and can thus potentially grow to amplitudes large 
enough to explain the spin equilibrium of LMXBs, but its 
modelling requires a detailed understanding of the physics 
of NS interiors. The r-mode can grow unstable if GW 
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emission drives it faster than viscosity damps it. This will 
only happen in a range of temperatures and spin frequencies 
which depends strongly on the details of the damping mech- 
anisms. In the standard picture, the main damping agent at 
low temperatures (below ~ 10 10 K) is the viscous boundary 
layer at the crust-core int e rface (|Bildsten fc Ushomirskvl 
|2000| ; iLevin fc Ushomirskvl |200ll l. while bulk viscosity 
is the strongest source of dampin g at high tempera- 
tures (|Andersson. fc Kokkotasl l200ll ). The nature of the 
damping mechanisms is very sensitive to the interior 
microphysics and presence of exotica, such as hyperons and 
deconfined quarks, or large scale superfluid and/or super- 
cond ucti ng components ( Andersson. Jon es & Kokkotas 
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r-mode oscillations distort the stellar magnetic field, lead 
to energy dissipation, and possibly prevent the mode f rom 
being driven unstable (|Rezzolla. Lamb fc Sh apiro 200(]). 

In this paper, we examine these mechanisms and com- 
pare them to observational constraints on NS spins and tem- 
peratures. We use available data on NS surface tempera- 
tures from X-ray observations of LMXBs in quiescence and 
UV observations of millisecond pulsars. We also present new 
analysis of five systems, which leads to new upper limits on 
their surface temperatures. We conclude that the minimal 
NS model, i.e., that of a star composed of neutrons, protons 
and electrons (possibly muons) and whose r-mode damping 
at low temperatures is due to Ekman pumping at the crust- 
core interface, is not consistent with observations and that 
additional damping mechanisms are required, unless the r- 
mode saturates at a very small amplitude. In this case GW 
emission would not affect the evolution of the system. 

We also discuss additional damping mechanisms that 
are likely to be at work in NS interiors and may be consistent 
with observations. 



2 R-MODE INSTABILITY WINDOW 

An r-mode is a fluid mode of oscillation of a NS for which 
the restoring force is the Coriolis force. To leading order in a 
slow rotation analysis, it is purely toroidal and has the form 



<5v 



I r V B 

a y— j RQY lm expiwt, 



(1) 



where <5v is the Eulerian perturbation of the total fluid 
velocity, Y^, is the magnetic-type vector spherical har- 
monic, R is the stellar radius and a is the (dimension- 
less) mode amplitude (|Owen et alj Il998l ). The fluid dis- 
placement gives rise to a current quadrupole moment and 
to the emission of GWs, which can drive the mode un- 
stable via the Chand r asekhar-Friedman-Schutz mechanism 
ilChandrasekharlll970l; I Friedman fc Schutzll 19781 ; lAnderssor] 



Il998l ; iFriedman fc Morsinkl Il998h . If GW emission drives 
the mode growth, then eventually the mode will saturate 
when energy is transferred to higher order modes due to 
non-linear couplings. Given the complexity of the full non- 
linear problem, this process is highly uncertain. Neverthe- 
less, most recent esti mates indicate a saturation amplitude 
10" 6 - 10" 5 , (jBondarescu. Teukolskv fc Wassermanl 



120071 ). This can be compared to an upper limit of a < 10 
from GW searches conducted with LIGO (|Owenll201(l ). Note 
that, as we are dealing with a superfluid star, the super- 
fluid neutrons can also flow independently from the charged 
component (protons and electrons), leading to relative mo- 
tion. In fact, to second order in the slow rotation analysis, 
the r-mode will ac quire poloidal components along the rela- 
tive velocity 5w„„ (IPassamonti. Andersson fc Has kell 2009; 
lHaskell. Andersson fc Passamontill2009| y 

An r-mode can be driven unstable as long as GW emis- 
sion drives the oscillation faster than viscosity damps it. 
This is usually studied in terms of the critical frequency at 
which the driving and damping timescales are equal. Solving 
for the roots of 



TGW 



TV 



(2) 



yields an instability curve that depends on frequency and 
temperature, tgw is the GW driving timescale which (for 
an I = m = 2 r-mode and an n = 1 polytrope) is given by 
jAndersson. fc Kokkot"asll200ll ) 



(3) 



with Mi. 4 is the NS mass in units of 1.4 M Q , i?io is the NS 
radius in units of 10 Km and P ms is the NS rotation period 
in milliseconds. The viscous damping timescale tv is given 
by 



(4) 



where the summation is over the various dissipative chan- 
nels, labelled with 'i'. At high temperature (above w 
10 10 K) the main contribution is bulk viscosity due to 
the modified Urea reaction, with a timescale given by 
l|Andersson. fc Kokkotasll200lf ) 



tbv = 2.7 x 10 11 M^aRT^pI.Tq 6 



(5) 



where T9 is the NS core temperature in units of 10 9 K. Note 
that this form for the bulk viscosity is only appropriate for 
small perturbations, such that perturbations of the chemical 
potentials is much smaller than the thermal energy kT. For 
much larger perturbations, the effect of bulk viscosity is sig- 
nifican tly stronger, effectively blocking the grow th of the r- 
mode l| Alford. Mahmoodifar fc Schwenzerl [20ll . However, 
the amplitudes that are necessary for such a scenario are 
significantly larger than the saturation amplitudes we con- 
sider here, so such a possibility will not be discussed further. 

At low temperatures, the main source of damping is 
the viscous boundary layer at the crust-core interface, which 
leads to a damping timescale 



TEK = 3 x 10 b P^fT 9 s, (6) 

where we use the estimate of [Glampcda kis fc Anderss"onl 
(2006) with a "slippage" factor S = 0.05. The slippage fac- 
tor accounts for the fact that the crust will not be com- 
pletely rigid, but will also participate in the oscillation. It is 
essentially the ratio between the crust/core velocity differ- 
ence and the mode velocity, so that S = 1 corresponds to 
a completely rigid crust, while smaller values indicate that 
the mode can penetrate the crust to some extent. Shear 
viscosity will also play a role at low temperatures, but its 
effect will be weaker than that of the crust-core interface 
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Figure 1. R-mode instability window for the "minimal" NS 
model described in the text, for which the main damping mech- 
anism at low temperature is the Ekman layer at the base of 
the crust. We schematically illustrate the trajectory a system 
would follow for high saturation amplitudes (a m 1) and low 
saturation amplitudes (a fa 10 -5 ). For large amplitudes, the sys- 
tem undergoes a thermal runaway and heats up significantly but 
spends much less than 1% of the time in the unstable region. 
For small saturation amplitudes, the time spent in the unsta- 
ble region increa ses, but the spin and temperature variation s are 
modest (fa 10%; [gondarcscu, Tcukolsky & Wasscrmanl l2007l ') . 



Furthermore if, as is generally believed, the core of 
the NS contains large scale superfluid components, these 
will rotate by forming an array of quantised vortices. 
The interaction of vortices with the charged components 
gives rise to a dissipative force known as mutual fric- 
tion. For temperatures well-below the superfluid transi- 
tion temperature, the damping timescale for mutual fric- 
tion is roughly constant, but the timescale can vary con- 
siderably when the temperature is near the transition 
temperature. Here we shall use the detailed results of 
lHaskell. Andersson fc PassamontH l|2009h . The main micro- 
physical input that is needed to calculate the mutual friction 
damping timescale is the value of the (dimensionless) drag 
param eter 71. It has been shown by Lindblom fc Mendelll 
(|2000h and lHaskell. Andersson fc Passamontil (120091 ) that 
the standard drag parameter (describing electron scatter- 
ing off vortex cores; 1Z ~ 10~ 4 ) does not significantly af- 
fect the instability window. However, the situation may be 
considerably different if the core of the NS is in a type II 
superconducting state. In this case, the magnetic field is 
arranged in flux tubes, and their interaction with neutron 
vortices could lead to strong dissipation, with drag param- 
eters possibly of the order of TZ sa 10~ 2 (| Jonesl fl992l ; iLinkl 
120051 ; lHaskell. Pizzochero fc Sidery||2oTll ). 

Finally it should be noted that, in magnetised stars, 
fluid motion distorts magnetic field lines, possibly leading to 
energy being drawn from the mode faster than GW emission 
can drive it (|Rezzolla. Lamb"fc Shapiro! [20001 ) . 



l| Andersson. fc Kokkotasll200ll ); thus we neglect it here. In 
Figure [TJ we show an example of an instability window for 
the "minimal" NS model described above and schematically 
illustrate the trajectory that an LMXB would follow. For 
large saturation amplitudes (a > 10 -3 ), the system under- 
goes a thermal runaway, i.e., heating up rapidly due to the 
unstable mode and t hen spinning down due to the emission 
of GWs (|Levinlll999l ). It would then enter the stable region, 
cool and spin up again, closing the cycle. For such a scenario, 
the system would spend less than 1% of the time in the un- 
stable regi on, making it very unlikely to observe a system in 
this stage (|Hevf 2002 V For more realistic values of the satu- 
ration amplitude (a ps 10~ 5 ), the system could spend up to 
30% of the time in the unstable region. The excursion in tem- 
perature and frequency are expected to be modest and one 
would still not expect to see a system well inside the insta - 
bility window l|Bondarescu. Teukolskv fc Wassermanll2007f ) . 
A further possibility is that the saturation amplitude is 
very small (which could be the case for low viscosity) and 
thus the GW torque is very weak. In this case, the GW 
torque cannot counteract the accretion torque, and a sys- 
tem would reach an equilibrium between r-mode heating 
and neutrin o cooling, while spinning up into the insta bil- 
ity window (jBondarescu. Teukolskv fc Wassermanll200^ . 

The minimal cooling model described above is modi- 
fied if there are hyperons in the NS core. The bulk viscos- 
ity can then be much stronger at lower temperatures. The 
effect of hyperon bulk viscosity on the r-mode instability 
window has been studied in detail for superfluid NSs by 
lHaskell fc Andersson! (|201Ch . whose results we shall use be- 
low. 



3 NEUTRON STAR TEMPERATURES AND 
SPIN RATES 

As is obvious from the discussion in Section [5] if we wish to 
construct an instability curve in the frequency versus tem- 
perature plane, it is necessary to estimate the temperature 
of the NS core, on which the damping timescales will de- 
pend. This is clearly not a straightforward task, as what 
is measured is the surface emission as detected by a dis- 
tant observer. In order to estimate the core temperature, we 
shall use X-ray observations of LMXBs in quiescence (when 
most of the thermal emission is thought to come directly 
from the NS surface; see, e.g.. lBrown. Bildsten fc Rutledgg 
1 1998T ) and the few available millisecond radio pulsar thermal 
spectra o bserved in UV. This is in c ontra st to the estimates 
made bv lHo. AnderssorTfc Haskell (2011), which made use 
of X-ray observations of LMXBs during bursts. 

Several LMXBs have surface temperatures obtained 
from blackbody fits to their observed X-ray spectrum. For 
others, the spectrum is completely non-thermal, and only 
upper limits on the temperature can be obtained. In table 
1, we list LMXBs that have a measured temperature (or up- 
per limit) and spin rate. The spin rates are either measured 
directly for those NSs that display coherent X-ray pulsations 
(indicated as "accretion powered") or inferred from the fre- 
quency of oscillations seen during thermonuclear type-I X- 
ray bursts (labelled as "nuclear po wered" ) . The s pin rates 
are taken from the overview given bv lPatrunol ll2010 | ). For the 
temperature s, we use the overview compiled by iHeinke al.l 
i|2007l, l2009h and include 10 additional sources reported in 
the literature or analyzed in this work (see Section 3.1). We 
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Table 1. Surface temperatures and spin rates for LMXBs that have measurements (or upper limits) of both. Note 
that Too is related to the surface temperature Tg by the relation T s = (1 + z)Too where 1 + z = (1 — 2GM / Rc 2 )~ 1 /' 2 
is the red-shift factor. We take 1 + z = 1.3 and the core temperature is then obtained as described in the text. 
We include three millisecond pulsars that have estimates of the surface temperature from UV observations. We 
distinguish between accretion powered (AP) pulsars (which are observed as X-ray pulsars), nuclear powered (NP) 
pulsars (for which the spin rate is estimated from type I X-ray burst oscillations) and radio pulsars (RPs). Pulsars 
that show both burst oscillations and X-ray puls ations are classified as accretion powered. Spin frequencies for the 
accreting systems are taken from iPatrunol feoiOT ). 



Source 


v (Hz) 


Too/10 6 K 


Tcore/10 8 K 


Type 


Reference 


Aql X-l 


550 


1.09 


1.08 


AP 


Heinke al. (2007) 


4U 1608-52 


620 


1.97 


4.55 


NP 


Heinke al. (2007) 


KS 1731-260 


526 


0.73 


0.42 


NP 


Cackett et al. (20101 


MXB 1659-298 


556 


0.63 


0.31 


NP 


Cackett et al. (20081 


SAX J1748. 9-2021 


442 


1.01 


0.89 


AP 


Heinke al. (20071 


IGR 00291+5934 


599 


0.82 


0.54 


AP 


Heinke al. (20091 


SAX J1808.4-3658 


401 


<0.35 


<0.11 


AP 


Heinke al. (20091 


XTE J1751-305 


435 


<0.82 


<0.54 


AP 


Heinke al. (20091 


XTE J0929-314 


185 


<0.58 


<0.26 


AP 


Heinke al. (2009) 


XTE J1807-294 


190 


<0.59 


<0.27 


AP 


Heinke al. (2009) 


XTE J1814-338 


314 


<0.80 


<0.51 


AP 


Heinke al. (20091 


EXO 0748-676 


552 


1.26 


1.58 


NP 


Degenaar et al. (2011) 


HETE J1900. 1-2455 


377 


<0.65 


<0.33 


AP 


This work 


IGR J17191-2821 


294 


<0.86 


<0.60 


NP 


This work 


IGR J17511-3057 


245 


<1.10 


<1.10 


AP 


This work 


SAX J1750.8-2900 


601 


1.72 


3.38 


NP 


Lowell et al. Qnjjreparation) 


NGC 6440 X-2 


205 


<0.37 


< 0.12 


AP 


This work 


SWIFT J1756-2508 


182 


<0.96 


< 0.78 


AP 


This work 


SWIFT J1749.4-2807 


518 


< 1.27 


< 1.61 


AP 


Degenaar et al. (in preparation) 


J0437-4715 


174 


0.12 


0.018 


RP 


Kargaltscv, Pavlov & Romani (2004) 


J2124-3358 


203 


<0.46 


<0.17 


RP 


Kargaltscv, Pavlov & Romani (2004) 


J0030+0451 


205 


<0.92 


<0.70 


RP 


Kargaltscv, Pavlov & Romani (2004) 



also include three millisecond radio pulsars for which the 
temperature was constrained by fitting the UV spectrum. 



3.1 Observations 

In order to obtain constraints on the surface temperature 
for some sources, we first simulate a fiducial X-r ay spec- 
trum with the software package XSpec [v 12.6; (|Arnaudl 
Il99rj )l. This is done using the NS atmosphere model NSAT- 
MOS (|Heinkeet al.ll2006r ), where we take M = 1.4 M 
and R — 10 km. We assume the entire NS surface is emit- 
ting (i.e., model normalization is fixed to 1) and use source 
distances reported in the literature (see below). After con- 
structing such a model for each source, we determine the NS 
temperature that produces the observed (quiescent) thermal 
flux limit. This value is then considered to be the upper limit 
on the NS surface temperature T s . 

We use flux upper limits reported in the literature to in- 
fer constraints on the surface tem perature for three source s: 
IGR J17191-2821 (D = 11 kpc; lAltamirano et aTll2010al l. 
NGC X-2 (D = 8.5 kp c; iHeinke et al l I 2010T I and Swift 
J1756-2508 (D = 8 kpc; iPatruno et alj|~201oh . In the case 
of IGR J17511-3057, nothing is reported in the literature 
about its quiescent properties. However, we found two obser- 
vations obtained with the X-ray Telescope (XRT) onboard 
Swift, which did not reveal the source during its quiescent 
state. We obtain an upper limit on the 0.5-10 keV unab- 
sorbed flux of ~ 7.5 x 10" 13 erg cm" 2 s" 1 . We infer a NS 
surface temperature of T s < 1.44 x 10 6 K for the fiducial 



mod el parameters mentioned above and assuming D — 6.9 
kpc l|Altamirano et a l. 2010b). We analysed a recent XMM- 
Newton observation of Swift J1749. 4-2807 during quiescence 
(Degenaar et al. in preparation). The source is clearly de- 
tected during the observation, but its X-ray spectrum is 
completely non-thermal . For a distance of D = 6.7 kpc 
ijAltamirano et al]|201ll l. we obtain an upper limit on the 
NS surface temperature of T s < 1.66 x 10 6 K. Finally, HETE 
J1900. 1-2455 has been continuously active since its discov- 
ery in 2005, but the source intensity dropp ed dramatically 
during a short ~ 20-day interval in 2007 l|Degenaar et al.l 
l2007d lbT). At a certain point the source could not be detected 
with Swift/XRT; t his resulted in an upper limit on the qui- 
escent X-ray flux (Deg enaar et al. 2007b). We re- analysed 
the data of this non-detection to estimate the upper limit 
on the NS surface tempe rature, where we assu med a source 
distance of D = 3.6 kpc (iGallowav et al.ll200Sl l. 



3.2 Neutron star core temperatures 

Having determined NS surface temperatures, we now esti- 
mate the core temperatures. We assume that the core and 
crust are nearly isothermal (which is very nearly the case 
since the thermal conductivity of the crust is high, as in- 
dicated by recent cooling observations of X-ray transients; 
I Brown fc Cummindlioog ) and that the core temperature is 
simply the temperature at the base of the heat blanketing 
envelope. As we are considering accreting systems and ra- 
dio pulsars that are thought to have been recycled through 
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accretion, we use the relation between surface temperature 
and enve lope base temperature for a parti a lly ac creted crust 
given bv IPotekhin. Chabrier fc YakovlevI (|l997t >. where we 
follow iBrown fc Cummind (|2009l ) bv "considering a layer 
of light elements down to a column depth of P/g = 10 9 
g/cm" 2 . 

We assume that millisecond radio pulsars have been re- 
cycled to rapid rotation by accretion and that they have 
similar crustal co mpositions as the LMXBs; thus we ag ain 
use the relation of lPotekhin. Chabrier fc YakovlevI (| 1997h to 
estimate their core temperature. This is of course a crude 
as sumption. However by using the i ron en velope relation 
of iGudmundsson. Pethick fc EpsteTnl |l983j), the estimated 
temperatures change by a factor of approximately 2. As 
we shall see, the MSRPs fall in a region of the instabil- 
ity window for which such a correction has no effect on 
our conclusions; this justifies our use of the relation from 
IPotekhin. Chabrier fc YakovlevI (|l997l ). Let us remark that 
the estimates of core temperature have large uncertain- 
ties, as not only is the composition of the envelope uncer- 
tain, but some systems may still be thermally relaxing to a 
steady state after an outburst and may have sizable temper- 
ature gradients in the crust. These effects lead to an uncer- 
tainty of a factor of a few in the inferred core temperatures 
(Bro wn fc Cummingl l2009V This has no qualitative impact 
on our conclusions, but this uncertainty should be kept in 
mind and we shall attempt to quantify it in the following 
sections. 



3.3 Strange stars 

It is possible that the most stable form of matter at the high 
densities that characterise a NS i nterior may be that of a 
conglomerate of deconfined quarks (|ltohll970l ; IWittenl 1984 ) . 
In fac t, it has been suggested that all NSs may be strange 
stars (|Alcock. Farhi fc Olintol Il986l 1. Although the ground 
state of matter at asymptotically high densities and low tem- 
peratures is known to be given by paired quar ks in the so- 
called "Colour Flavour Locked" (CFL) phase (|Alford et al.l 
2008), the properties of matter at realistic NS densities are 
still uncertain. The effect of a CFL core on the instabil- 
ity wi ndow was calculated bv lAndersson. Haskell fc Corned 
(2010) and found to be quite weak, so we shall not con- 
sider this possibility any further. We shall thus consider 
the effect of a conglomerate of unpaired quarks in the NS 
interior on the r-mode instability window, as this can be 
quite significant. In this scenario the shear viscosity dis- 
sipation timescale for a strange star (again assuming an 
n = 1 polytrope as the background model) is found to be 
l|Andersson. Jones fc Ko kkotas 20q3) 



(7) 



where as is the strong coupling constant. For the 
temperature range of interest (below ~ 10 9 K ) , the 
bulk viscosity damping timesca le is (Ma dsenl 1 19921; 
lAndersson. Jones fc Kokkotaslliool ^ 



tBv ~ 7.9 



/ m s y 
VlOOMeV/ 



Mi,4,R 10 T 9 

^ P-m s S 5 



(8) 



where m s is the mass of the strange quark. A strange star 
can also support a thin crust of normal nuclear matter 



up to the neutron drip density, after which free neutrons 
will drip into the strange core. Such a crust would obvi- 
ously be much thinner and much less massive than that of 
a NS, with a maximum mas s of approximately 10 Mq 
|Glendenning fc Webeij fl992). This crust would be much 
less rigid than a standard NS crust and not contribute 
signif icantly to the damping (|Andersson. Jones fc Kokkotasl 
120021 ) . We shall thus not consider damping due to the crust- 
core interface in the discussion on strange stars, although its 
inclusion would not qualitatively change our conclusions. 

The presence of the crust is, however, very significant 
for estimating the temperature of the core, as it provides a 
"heat blanket" for the strange core, allowing the outgoing 
radiation to thermalise (which would not be the case for a 
bare strange st ar, for which the sp ectrum would be consid- 
erably harder; iPage fc Uswl [2002'). Furthermore, many of 
the systems we consider show not only coherent pulsations 
but also thermonuclear bursts, which would be challenging 
to explain if there is no crust of normal matter. We shall 
thus use the same prescription as in the NS case to estimate 
the c ore temperature (see also the discussion in iPizzocherd 

mu). 



4 OBSERVATIONAL CONSTRAINTS ON THE 
INSTABILITY WINDOW 

Let us now examine how the observational evidence com- 
pares with theoretical calculations of the r-mode instability 
window. First of all, we begin by comparing the measured 
temperatures/spins with the minimal model instability win- 
dow of Section[2] In Figure[2]we show the inferred core tem- 
peratures and in Figure[3]we estimate the uncertainty due to 
the modelling of the outer layers of the star. The error bars 
on the core temperatures inferred from observations have 
been obtained by considering two extreme compositions for 
the stellar envelope, the properties of which (composition, 
thermal conductivity etc..) control the heat flow from the 
core to the exterior. We have thus calculated a "minimum" 
temperature by assuming a completely accreted cru st of 
light elements (|Potekhin. Chabrier fc YakovlevI Il997l ) and 
a "maximum" temperature by assumin g an iron envelope 
l|Gudmundsson. Pethick fc Epsteinl 19831 ). As we can see this 
produces an uncertainty of a factor of a few, which will dom- 
inate over the observational uncertainty but does not affect 
our conclusions. It is obvious from Figures [5] and [3] that, 
even accounting for theoretical and observational uncertain- 
ties associated with temperature measurements, several sys- 
tems are well inside the unstable region. As already men- 
tioned, this would is possible if the saturation amplitude 
is very large and one is lucky enough to catch the system 
while it is still in the GW emitting part of its duty cycle. 
However, given that for large amplitudes (a > 10 -3 ) one 
would expect the system to spend less than 1% of the time 
in the unstable region, it is highly unlikely that we are ob- 
serving so many systems in this phase. Furthermore, such a 
system would be spinning down rapidly due to the emission 
of GWs, but one of the systems, IGR J00291+5934, has a 
measure d spin-down rate in quiescence of 1) ~ 3 x 10~ 15 
Hz s - 1 llPatrunol 2010(: Hartman. Galloway fc Chakrabartvl 
120111 ; IPapitt o et al.|[201lf ), which is consistent with purely 
electromagnetic spin-down due to a B ~ 10 s G magnetic 
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Figure 2. R-mode instability window of LMXBs and MSRPs that have estimates of both the spin frequency and 
surface temperature (arrows indicate upper limits). The right panel is the same as the left panel but focused on 
the low temperature region in which the observed systems are located. It is obvious that a significant number of 
systems is well inside the "minimal" instability window, where one would not expect to find so many systems. In 
fact, for realistic values of the saturation amplitude, a star could not heat up enough to be significantly inside the 
unstable region, while for high values of the saturation amplitude a system would spend only a very small fraction 
of the time (less than 1%) above the instability curve, making it very unlikely to catch systems in this region. The 
only possibilities are thus that either the instability curve is significantly different from our minimal model curve 
due to additional damping mechanisms or the saturation amplitude is small enough not to affect the evolution of 
the systems. 



field (although one cannot rule out a much weaker mag- 
netic fi eld and low level GW emission, see lHaskell fc Patrund 
|201ll ) for a discussion of why this is unlikely to be 
the case in two other sources, SAX J1808. 4-3658 and 
XTE 1814-338). A final possibility is that systems in- 
side the window have undergone a thermal runaway and 
have reached an equilibrium between heating and cool- 
ing l|Bondarescu, Teukolskv fc Wassermanl [2003) ; they are 
now either at spin equilibrium (i.e., with the GW spin- 
down torque balancing the accretion spin-up torque) or ap- 
proaching spin equilibrium (as could be the case for IGR 
J00291+5934 which exhibits long-term spin- up). We discuss 
this possibility further in the following section. It is, how- 
ever, clear that the minimal model is not consistent with 
observations. 

We now discuss the possible mechanisms that may be 
at work in a realistic NS and that could be consistent with 
observations. We first examine effects due to properties of 
the crust. One is that the crust may be more rigid than 
is commonly assumed. This would lead to stronger dissipa- 
tion at the crust-core interface. In Figure 3J we show the 
effect of increasing the "slippage" factor S from a stan - 
dard value of S — 0.05 (|Glampedakis fc Anderssoni 2006) 
to S = 1 (a completely rigid crust). It is obvious that a 
mor e rigid crust could allow all the systems to be stable (see 
also I Wen et a.1.1 T201 lh . However such a rigid (<S = 1) crust 
is not realistic. Alternatively, given the frequency range of 
r-modes, the mode may couple effectively to torsional os- 
cillations of the crust. This wo uld produce strong dissipa- 
tion at the resonance frequency (iLevin fe Ush omirskv 200 ll; 
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Figure 3. The same r-mode instability window as in the right 
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due to the uncertainty in modelling the outer layers of the NS, 
as described in the text. We can see that although there is a 
significant uncertainty on the inferred core temperatures it is not 
large enough to modify the conclusion that many of the systems 
appear to be well inside the unstable region. 



Another possibility is that core bulk viscosity may be 
much stronger at low temperatures. For example, if hyperons 
are present in the core, then a significantly restricted unsta- 
ble region is created, as illustrated in Figure [5] The situation 
for strange stars is somewhat similar, with bulk viscosity 
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Figure 6. R-mode instability window for a strange star with 
M = 1.4M0 and R = 10 km. We assume unpaired quarks in 
the core (see text). The observed temperature would require the 
strong coupling constant to be smaller than generally assumed 
and a rather large strange quark mass. 
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Temperature (K) 

Figure 5. R-mode instability window for a M = 1.4A/0 
star with hyperons in t he core. We use the results of 
lHaskell fc AnderssorJ l|2010h for different NS radii and values of 
the coupling parameter v , whic h parametrises in- medium effects 
(sec Haskell Sz Andcrsson (2010) for details). Note that the inter- 
ruption at low temperatures for the \ = 1 curve is not physical, 
but due to difficulties with the numerical setup. It would appear 
that a rather compact star and strong coupling are required for 
the model to be consistent with observations. 



playing a much stronger role at low temperatures and lead- 
ing to a reduced unstable region, as show n in Figure [6] (see 
lAlford. Mahmoodifar fc Schwenzerlkoill . for more detailed 
discussion of the r-mode instability window in strange stars 
and hybrid stars). Although interesting, this mechanism has 
several pro blems. First, the meas urement of a NS with mass 
hi ~ 2Mq l|Demorest et alj|2010l) can already exclude so me 
hyperonic equations of state ( Lattimer fc Pr akash 201(3; al- 
though the presence of hyperons in the core may be consis- 



tent with the low observed te mperatures of som e sources re- 
ported in lHeinke al"1l2007l and lHeinke"al1l2009h . Second, one 
cannot explain the existence of millisecond radio pulsars: af- 
ter accretion has ceased, one would expect most systems to 
cool and pass through the unstable region, during which the 
NS would spin down and not maintain the high spin rates 
that are observed (y w 600 — 700 Hz). Finally one needs hy- 
peron/quark bulk viscosity parameters that are somewhat 
extreme to reconcile with the observed temperatures and 
spins. 

A very interesting mechanism is one that involves strong 
vortex-mediated mutual friction. If the core of the NS con- 
tains a type I superconductor, then mutual friction will not 
be st rong enough to significantly affect t he in stability win- 
dow (|Sedrakianll2005l ; although see |jonesil2006l . for a discus- 
sion of strong drag in type I superconductors). However, if 
the core contains a type II superconductor, then the inter- 
action of vortices with flux tubes will lead to strong mutual 
friction if a large fraction of vortices can "creep" through 
the flux tubes. Examples of this are shown in Figure [7] 

Finally a pro mising scenario involves magnetic damp- 
ing of the r-mode (|Rezzolla. Lamb fc Shapiroll2000l ). Given 
the high electrical conductivity of the NS interior, the mag- 
netic field lines are frozen in with the fluid and can thus 
be distorted and wound up by the oscillatory motion of the 
r-mode. Even for relatively weak magnetic fields, this could 
lead to rapid damping an d could close the instabi lity window 
IjRezzolla et al.ll2001al lbt ICuofano fc Dr ago 2010). 



5 SPIN EQUILIBRIUM 

We now examine the possibility that GW emission due to 
an unstable r-mode may be setting the spin equilibrium for 
LMXBs. This could be the case if the critical frequency 
increases with temperature at around 10 7 K (e.g., for hy- 
peron and quark bulk viscosity or for strong mutual fric- 
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tion). As a result, thermal runaway is halted, and the system 
reaches an equilibrium state, such that viscous heating due 
to the r-mode is balanced by neutrino emission and the GW 
torqu e balances the accretion torque at the observed spin pe- 
riod l|Andersson. Jones fc Kok kotas 20021; iNavvar fe Owen] 

l2ooeh . 

W e follow the approach of I Ho. Andersson fe Haskell 
and assume that a GW torque due to an unstable 
r-mode is balancing a spin-up torque due to accretion on 
the N S surface. The heat dissip ated in this case has the 
form (|Brown fc Ushomirskv||2000r i 



L heat= 0-064 (^-) L a 



(9) 



Taking the heat from r-mode dissipation to be lost by neu- 
trino emission [i.e., i nea ^ = L V (T), where L v is the neu- 
trino luminosity], the core temperature T can be inferred. 
In order to determine the rate at which neutrino emission 
cools the system, it is important to account for superfluid- 
ity, as this will lead not only to a reduction in the emis- 
sion rates for the modified Urea emission processes but 
also to additional neutrino emission from the formation of 
Cooper pairs. We use the latest constrain t s on superfluid 
trans ition temperatures (|Page et all l201ll ; IShternin et al.l 
l201ll ). obtained from the observed rapid c ooling of the 
NS in the Cassiopea A su pernova remnant (|Heinke fc Hoi 
2010: Shtc rnin et al.l |2011|). We use the code described in 
Ho. Glampedakis fc Andersson! (|201lT ) to calculate the neu- 



trino luminosity, which is obtained by integrating the neu- 
trino emissivities over the stellar volume. Briefly, this in- 
cludes building a NS with the APR I equation of state (in 
our case, M = 1.4 Mq and R = 12 km), calculating neutrino 
emissivities due to modified Urea, nucleon bremsstrahlung, 
Cooper pair formation and breaking, plasmon decay, and 
pair annihilation in the core and crust, and accounting for 
neutron singlet and triplet superfluids in the core and crust, 
respectively, and proton superconductivity in the core. Note 



that th e difference between core temper atures derived here 
and in I Ho. Andersson fc Haskelll ()201 lh is < five percent 
and the difference between core temperatures derived us- 
ing ne utron triplet transition tem peratures from lPage et al.l 
J201ll) and Shternin et all J201ll ) is < fifteen percent; see 
iHo. Andersson fc Haskelll l | 20 lift , for derived core tempera- 
tures assuming only modified Urea neutrino emission. 

In Figure O we show the temperatures obtained in the 
spin equilibriu m scenario for a "shallow" neutron su perfluid 
transition fsee lHo. Glampedakis fc Andersson ( 2011 ) for de- 
tails) with Tm m „ sa 5 x 10 s K. as in lPage et all (l201ll). The 
long-t erm accret ion luminositie s are t aken from I Watts et afl 
l|200Sft and from lFalanga et al l (|201lft for IGR J17511-3057. 
We can see that many systems appear to be colder than what 
would be expected in the presence of an unstable r-mode, 
although for some of the faster systems (which are also the 
most likely targets for GW searches, given the strong scaling 
with frequency of the GW torque), GW-driven spin equilib- 
rium may still be possible and cannot be completely ruled 
out. 

Finally we can calculate the maximum amplitude that 
would be compatible with the inferred core temperatures by 
assuming that the viscous heating is due to an unstable r- 
mode with arbitrary amplitude a ( jAndersson. fc Kokkotas! 
l200lft 



^heat 



1.31 



a 2 v 2 MR 2 



(10) 



where M is the mass of the star, R its radius, we have as- 
sumed that the equation of state is given by an n — 1 poly- 
trope and r av is the shear viscosity damping timescale. If 
we assume that shear viscosity is mainly due to electron- 
electron scattering and that modified Urea reactions are the 
main contribution to the cooling, the maximum amplitude 
we obtain takes the form 

n5 



6.7 x 10" 5 — 
v 



(11) 
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Figure 8. The left hand panel shows the core temperatures derived from X-ray luminosities of LMXBs measured 
during bursts and assuming spin equilibrium (diamonds), and core temperatures derived here from UV and quiescent 
X-ray luminosities (triangles). The panel on the right hand side is the same but also shows the range in temperatures 
obtained by assuming that the minimum core temperature is given by assuming an envelope completely composed 
of light elements, while the maximum is for a pure iron envelope, as described in the text. It is obvious that many 
systems are too cold to allow for a spin equilibrium r-mode. However, the more rapidly rotating systems, which are 
hotter, may be consistent with spin equilibrium. 



where T% is the temperature in units of 10 s K and v is the 
spin frequency of the system. With the notable exception of 
the two fastest and hottest systems, 4U 1608-52 and SAX 
J1750. 8-2900, these are more stringent upper limits on the 
mode amplitude than the spin equilibrium condition and 
from the results in table [1] we see that for most systems in 
the unstable region one has a m « 10 -9 — 10 -8 , which could 
be compatible with a small saturation amplitude for the r- 
mode a nd would have no impact on the spin evolutio n of the 
system (jBondarescu. Teukolskv fc Wassermanll2007l) . 



6 CONCLUSIONS 

In this paper, we estimated the core temperature of NSs 
using data from X-ray observations of LMXBs in quiescence 
and UV observations of millisecond pulsars, in order to place 
constraints on the physics of the r-mode instability window. 
We also presented a new analysis of five systems. 

These estimates show that, if one uses a "minimal" NS 
model, in which shear viscosity is due to dissipation in a 
boundary layer between the crust and core and bulk viscos- 
ity is due to modified Urea processes, the r-mode would be 
unstable in many of these systems. In particular, many sys- 
tems are well above the critical frequency-temperature for 
the instability to grow, which is highly unlikely since systems 
should never depart significantly from the stable region for 
sm all saturation mode amplitudes (such as tho se predicted 
by iBondarescu. Teukolskv fc Wassermanl 120071 ) , while sys- 
tems should spin down too quickly to be detected for much 
larger values of the saturation amplitude (Heyl 2002). 

It is clear that additional physics and additional damp- 
ing mechanisms have to be built into the model for it to be 
consistent with observations. Enhanced bulk viscosity, due 
to hyperons or deconfined quarks in the core, could provide 



a source of damping that is consistent with observations, but 
their presence would also predict that, as systems cool after 
accretion ceases, they should once again enter an unstable 
region and spin down. This is at odds with the existence 
of rapidly rotating millisecond radio pulsars. Strong mutual 
friction may also be consistent with observations if super- 
fluid vortices can cut through superconducting flux tubes 
(see Glampedakis, Andersson & Haskell, in preparation, for 
a detailed discussion of this scenario). 

A promising scenario is that in which damping is due 
to the crust responding rigidly to the mode displacement in 
the r-mode frequency rang e for rapidly rotating system. The 
phenomenological model of lHo. Andersson fc Haskelll (|201ll ) 
shows that this is viable, but more quantitative models are 
needed. In particular, efforts should be made to better un- 
derstand the effect on viscous da mping timescales of pasta 
phases at the crust-core interface l|Horowitz fc Berrv||2008D 
. Another scenario that would be consistent with observa- 
tions is that in which the mode winds up the magnetic 
field of the star, and energy is extracted from the oscilla- 
tory motion more rapidly than GW emission can drive it 
l|Rezzolla. Lamb fc Shapiroll2000l) . Once again this mecha- 
nism depends strongly on the internal magnetic field struc- 
ture and further work is needed in order to assess its rele- 
vance, as well as accounting for the presence of supercon- 
ducting components. 

Finally an interesting possibility is that the satura- 
tion amplitude of the r-mode is small enough that the 
GW torque cannot counteract the accretion torque and 
a system would spin up into the unstable region. In or- 
der for this scenario to be consistent with observations 
(i.e. in order for the heating from the mode to be consis- 
tent with the observed temperature), the saturation am- 
plitude should be roug hly a < 10 -9 - 10~ 8 . Such a small 
amplitude may be consistent with theoretical calculations 
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l|Bondarescu. Teukolskv fc Wassermanll2007l ') and would in- 
deed lead to a spin-down torque that is smaller than the 
electromagnetic spin-down torque for a B ~ 10 8 G magnetic 
field, thus not impacting on the evolution of the systems. 

We examined the possibility that continuous GW emis- 
sion from an unstable r-mode may be setting the spin equi- 
lib rium period of the LMXBs. T his scenario was considered 
bv lBrown fc Ushomirskvl |2000l ) who found that, if one as- 
sumes modified Urea cooling, most systems would be too 
hot to be consistent with observations. We re-examined this 
scenario by using the most recent constraints on superfluid 
transition temperatures obtained from observations of the 
cooling of the NS in Cassiopeia A. We find that this leads to 
lower core temperatures (due to stronger neutrino emission) 
which may be consistent with the more rapidly rotating sys- 
tems. This is interesting since the GW spin-down torque 
scales strongly with frequency and is expected to play a 
stronger role in rapidly rotating systems. Further observa- 
tional constraints, as may be available from future X-ray 
observatories such as LOFT and Astrosat, as well as theo- 
retical work on NS composition and viscosity, are crucial to 
aid in the search for GWs from these systems (|Watts et al.l 
l200Sh . 
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